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Abstract 



Tlie dynamics of an idealized, infinite, MIT-type flux tube is followed in 
time as the interior evolves from a pure gluon field to aq q plasma. We work 
in color U(l). q q pair formation is evaluated according to the Schwinger 
mechanism using the results of Brink and Pavel. The motion of the quarks 
toward the tube endcaps is calculated by a Boltzmann equation including col- 
lisions. The tube undergoes damped radial oscillations until the electric field 
settles down to zero. The electric field stabilizes the tube against pinch insta- 
bilities; when the field vanishes, the tube disintegrates into mesons. There is 
only one free parameter in the problem, namely the initial flux tube radius, 
to which the results are very sensitive. Among various quantities calculated 
is the mean energy of the emitted pions. 
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I. INTRODUCTION 



Flux tubes are one of the most elementary systems of quantum chromodynamics. They 
are the idealized configurations of heavy quark-antiquark pairs at large separation L such 
that the region between can be assumed to possess axial-cylindrical symmetry. They play a 
central role in lattice QCD calculations and in models of QCD, as well as in the phenomenol- 
ogy of QCD processes, such as jets produced in heavy ion collisions or the "spaghetti" 
connecting the separating heavy ions. 

From the spectroscopy of heavy quarkonia and Regge trajectories, it has been possible 
to extract a reliable measurement of the string tension by assuming the two body potential 
for heavy quarks to be of the form 

V{L) = ^ + 91+ constant (1) 

The "experimental" value is Jll ^ = 913 MeV/fm = 0.18 GeV^=4:.63/fm^. In fact the linear 
region of the flux tube potential is not directly explored by quarkonia or by lattice QCD 
calculations: 

The wave functions in quarkonia calculations tend to be centered near the "knee" of the 
potential although high-lying states do explore more deeply into the linear region. However 
other forms of the heavy quark potential (including a logarithmic term) have also achieved 
some success in reproducing quarkonia spectra. 

Lattice QCD calculations on flux tubes are generally limited to the quenched approxi- 
mation (no massless quarks) and allow for a separation of the heavy quark-antiquark of only 
about 1 fm, which do not reproduce a cylindrical region 

Static flux tubes are unstable at separations greater than 1 fm, since the energy required 
to stretch the tube by 1 fm is nearly 1 GeV, which is about the energy difference between 
a quarkonium, QQ, and a pair of heavy-light mesons, Q q + Q Q- Lattice calculations 
without light quarks cannot explore this instability. 

Recently there have been several papers which have explored the creation of light 

quark-antiquark pairs as a mechanism for flux tube breaking. They employ a variation of 
the Schwinger parallel plate capacitor model in which the infinite transverse geometry is 
replaced by the MIT boundary condition 
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(2) 



at the cylindrical tube surface for the light quarks. The electric field E is constant over all 
space in the longitudinal {z) direction, with D equal to E inside the tube and zero outside. It 
is found that pair production in the cylinder is suppressed relative to the Schwinger formula 
due to transverse confinement which gives the quarks an effective mass proportional to 1 /a 
where a is the radius of the flux tube. 

In the present work, we make some simplifying assumptions, the most significant are: 
We consider infinite flux tubes. This means that the tube is long compared with its 
radius. For a jet, this means that the characteristic tube lifetime, r >> a/c, since the quark 
producing the jet travels near the speed of light. 
We use the MIT bag model for the tubes. 

Although some contact is made with color SU(3), we perform our calculations in U(l). 
The extension to Abelian SU(3) is not a serious complication. 

We adopt one dimensional Debye screening for the quark-quark interactions used in the 
Boltzmann collision term. This is further simplified to an effective (5-function potential. 

With these approximations, our results depend on only one parameter, the initial (static) 
flux tube radius. 



II. THE MIT MODEL 



We utilize the framework of the MIT model because it is so very simple. It is the proto- 
type of color dielectric models. Some features of various models of flux tube are discussed 



in Sec. [VITI . 

In the MIT model, the calculation of the flux tube configuration is straightforward since 
it involves no quarks. The flux through the tube is equal to the color charge Q = ^\gs, 
where the A are the Gell-mann color matrices. The electric field is E = Q/A. Here Q is the 
charge at one end (say, the left) of the flux tube {Q at the other). It is this charge which 
will be shielded as a function of time due to the uniform production of pairs out of the sea 
in the flux tube. We distinguish it from the charge q of the light quarks being created in 
pairs, although for the static model Q = q. 
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The energy per unit length is given by 

E 1 

where A = vra^, with a the cyhnder radius, and B is the bag constant. Minimization with 
respect to A yields A = Q/ \J2B and 



= e 



2 



2Q'b]'^' = 2BA = ^. (4) 



L 

The translation to QCD obtains by the replacement of by its expectation value in a 
color-singlet state, 

0^ >= I < >= M^, (5) 
where we have used < >= 16/3 and as = fi's/47r. 

III. THE SCHWINGER MECHANISM 

The rate of pair creation in the Schwinger mechanism with MIT boundary conditions 
has been calculated by several groups 0, 0, [Q. The point of departure is the QED 
Schwinger problem of infinite parallel plate capacitors infinitely far apart 0. The rate of 
pair production per unit time per unit volume is 



^2 oo ^—Trnm'^/K 
n=l 



for mo — > 0. (6) 



47r3 6 247r 

where mo is the Fermion mass and k = IgEl. 

In the MIT model, k (with q = Q = VOna) is just the string tension 6. in which case 

w{oo) ^ 0.25 c/fm^ (7) 

This identification is model-dependent, as discussed below. For example, Pavel and Brink 
identify n = 26. (For reference, we note that those authors denote the string tension by 
cr whereas we use 9.) 

As we will see below, the rate is strongly suppressed by confinement. 



IV. CONFINEMENT AND TRANSVERSE MASS 



Transverse confinement has the effect of introducing a "transverse mass" in the formula 
for w. rriQ itLq + x^/a^. Using the results of Pavel and Brink [cf. their Eq. (4.3)], the 
pair production rate per unit time per unit length W is 



W{a) 



71 



Xn>0 



1 



-TTX^/a^K 



nrK 



-■Kxl/a'^K 



TT 



(8) 



where is the number of flavors. For reasonable a and itlq = 0, the two lowest transverse 
modes, x±i dominate, and we consider these only here: \x±i\ = 1.4347. Henceforth we will 
assume uiq = 0. 

V. SIMPLE ESTIMATE OF COLLAPSE TIME 



In the following, we work in QED, which is to say, charge is a scalar quantity, plus or 
minus. 

We now consider Q = Q{t), since the charge at the ends is shielded due to the flow of 
opposite charge into the ends. We will calculate it more carefully in the next section, but it 
is a good approximation to calculate the current by using the initial radius of the flux tube 
and the electric field for short times. The number of pairs is proportional to the time, and 
the mean velocity of a quark is proportional to the time. It is then easy to show that 



-2 J 



-2 



(9) 



dt m 2 Am 

where J is the current carried by positive quarks; negative quarks carry an equal current, 
including the same sign. Hence the factor of 2. We have used the average velocity of quarks 
< p > /m ^ Q{0) qt/{2Am), We make the identification m = m± = xi/a. Again we note 
that a is the tube radius at t = 0. One more time integration gives 



Q{t) = g(o) 



3 Am 

We now estimate the characteristic lifetime of the tube to be 



(10) 



ria] 



3mA 


1/3 


'3m' 


1/3 


"3 7rxe''^i/'*''^ 


_Wq\ 




[wk\ 




a 



1/3 



(11) 
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Note that there is no dependence on the length of the flux tube! 
Taking = 2, we find 



r(0.5) = 7.0 
t(1.0) = 1.4 

Recall that k = 9 in the MIT model. 



29 (12) 
1.7 fm/c, (13) 
0.68 fm/c. (14) 



VI. COLLISIONLESS BOLTZMANN EQUATION 

Here we still employ the adiabatic MIT model to follow the dynamics of the flux tube 
wall by assuming that the tube radius is in instantaneous equilibrium. In addition to the 
electric field, the light quark plasma also exerts a pressure on the tube wall. We assume that 
only the lowest transverse eigenmode is excited {x±i), so that the tube radius, now denoted 
by R{t) instead of a, is determined by minimizing the transverse energy 

Z = "^^+2^ + ^' 

which yields 



TT 



2 



+ 27rBR^ -2NxiR = Q. (16) 



N — N{t) is the number of pairs per unit length, 

n{p,t)dp= / W{t')dt (17) 

-OO Jo 

with n(p, t) the momentum density per unit length for either qorq. Note that we distinguish 
between the shielded heavy quark charge Q{t) and the elementary light quark charge q, where 

g(o) = q. 

The pair production rate (per unit length) is given by 

= -6-'^"'/^'% (18) 

TT 

with K, = q\E\. 
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The transverse motion of the quarks is now treated relativistically but classically. Thus 
the relativistic mass is taken to be 



Xi 



2 



'"^("^jgyj -^^ ■ (19) 

The rate of change of n{p,t) is governed by the coupled integro-differential equations: 

Ji^) = Q / dpnp/m, (21) 



§^-2 J, (22) 



where {dn/dt)coi is the rate of change of the momentum distribution due to collisions. The 
initial conditions are n(p, 0) = 0, Q{0) = q = VOttu, -R(O) = a. 

We wish to point out the very interesting fact that no magnetic fields are produced! 
This is because, in Maxwell's famous fourth equation, the convection current is exactly 
cancelled by the displacement current. 

One valuable check on the numerics is provided by the agreement between the two forms 
of Eq. Excellent agreement is attained until numerical instability sets in, when the 
two values diverge dramatically, in which case time or momentum steps were changed until 
equality was achieved. 

These equations have been solved numerically first without the collision term for several 
parameter sets. The results are displayed in figure ^ the thin lines are for the collisionless 
solutions. We note here oscillations with negligible damping. In the absence of damping, the 
current continues to flow even when the charge at the end caps (and hence E goes through 
zero. Consequently the charge at the end caps changes sign, tending to reverse the current. 
The number of pairs increases with time. The energy per unit length also increases; there 
is no energy conservation within a segment of the tube. 
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VII. COLLISIONAL THERMALIZATION AND DAMPING 



The problem is simplified by the one-dimensional nature of the dynamics: we assume 
that the transverse motion can be handled by the adiabatic MIT model. Scattering between 
quarks (charge q) or between quarks and anti-quarks (charge —q) is like that of beads on a 
string. The colliding particles exchange momentum with some probability to be computed 
below. Collisions between identical particles does not alter the momentum distribution - 
it is as though no collision took place. Since the interaction is independent of spin and 
flavor, we do not consider any quark-quark collisions, but only collisions between quarks 
and antiquarks. Then the change in population density of quarks in momentum satisfies the 
equation 



with n = n_|_. Here |7^p is the reflection probability for the collision of quarks with antiquarks 
in a one dimensional potential barrier or well. The potential should take into account 
shielding introduced by the presence of other quarks and antiquarks. 



The interaction between confined quarks is that of disks in the tube of radius R. The 
confinement of the D-field by the tube walls leads to a potential between bare quarks which 
is linear in the separation, just like the flux tube potential itself. However the quarks are 
shielded by other quarks and antiquarks. An estimate of the shielding can be obtained from 
a one-dimensional version of the Debye formula: 

The charge density in the vicinity of a quark at (say) 2; = in a thermal bath in the rest 
frame of the quark is 




dp2 |7^(pi -p2)f\{vi - V2\[n-{pi)n+{p2) - n+(pi)n_(p2)] 




(24) 



A. The Shielded Quark- Antiquark Interaction 




(25) 



where 




(26) 
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with mo = Xi/R and (pi^z) the shielded potential to br determined. Since 0(z = ±00) = 0, 
it follows that the polarization charge of +/— quarks is 



p^{x = ±00) 



and hence 



q .( ^ , qN 



±C I e-V^^/f^BT^p 
J —00 



(27) 



(28) 



If, as turns out to be the case generally, \q(t))\/kBT << 1, the expression can be linearized 
in 0, and we obtain 

2Nq'^ 



q 



■5{z) 



Poisson's equation in one dimension then reads 

= T^^i^) + 



The solution to this equation is given by 



nR^kBT 



q\ 



with the Debye length given by 

2 _ nR^keT 
2iVg2 

Therefore the quark-antiquark potential is given by 



Vn(z) 



-q(j){z) 



q^^D \z\/Xo 

27rR^ 



(29) 



(30) 



(31) 



(32) 



(33) 



B. Calculation of the reflection probability 

It is an elementary quantum mechanical problem to calculate |7?-|2 for the one dimen- 
sional Debye potential Eq. For greater simplicity and insight, we consider two analytic 
approximations for |7?.|2: 

(1) A square-well potential of width D and depth Vq such that VqD equals the integral 
of Vd{z). a reasonable value for D is 2A/). The reflection probability can be found in any 
of several texts. 
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(2) A delta-function potential of strength Vs (like the square- well potential with D 0). 
The strength is given by 



C. Evaluating the temperature 

In order to identify the temperature T, we note that, except for very small times, the 
charactersitic momenta in the problem are large compared with m, and therefore we may 
set ~ ±q(j) + \p\ in Eq- ^ Consequently we can approximate [cf. Eq. |25[ 

n{p) = n+ip) ^ Ce-IP~<*'>l/'^^^, (35) 

where < p > is the average drift momentum. For this form of the Boltzmann distribution, 
we have 

/oo 
n{p)dp = 2CT (36) 
-oo 

and the fluctuation in the momentum yields the temperature: 

/oo 
n{p)\p- <p>fdp (37) 
-oo 

In the calculation, we use these equations to determine the temperature at each time 
from the calculated n{p,t). Then 

{kBTf =< {p- <p>f> /2 (38) 

where < > means momentum average with respect to the calculated n{p,t). 

D. TZ for the delta- function 

The delta-function case V{zi — Z2) = VsS{zi — Z2) has a particularly simple form for the 
reflection amplitude TZ and probability |7^p for particles of momenta pi and p2, namely 

^ W+Jp^^WJ? ^^^^ 
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where, without rigor, we have used a reduced mass 

111 

fi m{pij m[p2) 

With the self-consistent determination of Vs from Eqs. ^ and |^, a complete solution of 
Eq. ^ and comparison with the collisionless results for the same model parameter (only a) 
is possible. 



E. Damping time 

We see here the evolution from a flux tube to a plasma tube. The behaviour of this 
evolution depends strongly on the single parameter a. For a between 0.5 and 1.5 fm, we 
note damped oscillations which increase in amplitude and decreases in period as a increases. 
We proffer three measure of the characteristic damping time, tq, ti, and T2. tq is defined as 
the time at which Q{t) first crosses zero. The other definitions are 

n = / \Q{t)\tdt/ \Q{t)\dt (41) 

"'0 / 

and 

T2 = J^ Q\t)tdtl Q{tfdt (42) 

Calculated results for these three measures are shown in Fig. |[ 

In Fig. ^ we show the number of pairs produced per unit length, and the energy per 
pion. The pion production is calculated assuming that one meson, vr, p, or uj, is produced in 
proportion to their spin-isospin statistical weights: 3/15 for pions, 9/15 for rhos, and 3/15 
for omegas. The rho and omega decay to 2 and three pions respectively. The total energy 
per unit length is distributed among the pions per unit length. The result is shown in Fig. 
i 

Note that our calculations do not depend upon the motion of the end caps, and therefore, 
within our approximations, the energy per ion is flat in rapidity. 

Our results show a high degree of sensitivity to the initial flux tube radius a. This 
provides a useful handle for distinguishing among various flux tube models and, in particular, 
the appropriate size of a. 
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F. Ultimate disintegration of the plasma tube 



Unlike the plasmas in a fusion reactor, the flux tube is stable against the sausage insta- 
bility; the difference is due to the absence of magnetic fields. After the electric field had 
decayed, this stabilizing effect is removed and the plasma tube falls apart. We have not 
attempted yet to follow this final phase, but use the various r's given above to estimate the 
disintegration time. Our own bias is to use ti. 

VIII. MODELS OF FLUX TUBES 

Since the collapse time is strongly dependent upon the initial flux tube radius -R(O) = a, 
we now address more realistic flux tube calculations to determine a and the relationship 
between n and 9. 

A. The MIT Model 

In this model, the field energy is uniform within a radius a which is characteristically 
1 fm or larger. The flux tube energy is shared equally between the electric field and volume 
energies: 9 = E^A = Q^/A = q\E\. 

B. Lattice Gauge Calculations 

Lattice calculations at present are usually restricted to the quenched approximation (no 
zero-mass dynamical quarks) and to Q Q separations L of less than about 1 fm. Nevertheless 
Sommer 0] obtains a dependence of the (longitudinal) electric field energy which appears 
to be stable with L. The functional form is roughly exponential, E'^/2 oc e^'^'i/'' with 
< r\ >= 6 6^ ~ (0.2 fm)^. The equivalent radius (the radius of a square form giving the 
same < rl_ >) is Veq = V2 < >^/'^^ 0.3 fm. Here 9 = 2q \E\. 
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C. Dual Superconductivity Model 

The dual superconductivity model yields a somewhat larger radius than lattice cal- 
culations, with Ten =~ 0.4 fm. 



D. The CDM 



The chromodielectric model is similar to the MIT model, except that the electric field 
energy has a diffuse surface and the flux tube energy 6 also contains a surface term such 
that e > q\E\. 



IX. FURTHER CORRECTIONS 



There are several effects which are yet to be included in order to improve the estimate 
of the collapse time: 

(1) Finite length effects also increase the collapse time. There are at least two such effects 
to consider: a) The current associated with a given momentum group terminates when all 
the charge in the group is transported the full length of the tube, b) No pairs can be created 
unless the change in potential the length of the tube is greater than m^: \qE L\ > Xi/R. 
These also lengthens the collapse time. 

(2) Dynamics of the conflnement mechanism and associated energy conservation. One 
might consider the chromodielectric soliton model to handle this. 



X. FINAL COMMENTS 



It is important to note that the calculations presented here do not depend on the velocity 
of the heavy quarks in the longitudinal direction. The collapse time depends sensitively on 
the radius of the flux tube and indeed on detailed structure of the tube. 

For different approaches to the problem, see Kluger, et al. [|rU| and Rau [O . 
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FIGURES 

FIG. 1. Schematic evolution of a flux tube to a plasma tube to a collection of mesons. The 
radius of the tube exhibits damped oscillations. 

FIG. 2. Variaton with time of the flux tube radius R{t)/a (upper curves) and net 
endcap charge Q{t)/q (lower curves) for various a. For the case (a) of a = 0.5 fm. 
we have also shown in thin lines the solution to the collisionless Botzmann equation. 

FIG. 3. The mean longitudinal momentum per quark as a function of time for a = 0.5 and 1.0 

fm. 

FIG. 4. The quark "temperature" as a function of time for a = 0.5 and 1.0 fm. 
FIG. 5. The growth of the number of pairs as a function of time for a = 0.5. 
FIG. 6. The variation of several important quantities with a: The characteristic damping times 



(in fm/c) To (where Q{t) flrst crosses zero), ri and T2 (as defined by equations 41 and |4^) . Ten 
times the number of pairs produced per fm is denoted by 10 N. E is the final energy per fm (1/fm^). 
£^7r is the mean energy (in GeV) of pions produced by disintegration of the tube. 
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